arXiv:l507.00160v2 [cond-mat.str-el] 7 Jul 2015 


EPJ manuscript No. 

(will be inserted by the editor) 


Extended s-wave pairing symmetry on the triangular lattice 
heavy fermion system 

Lan Zhang^^, Yu-Feng Wang^, Yin Zhong^, and Hong-Gang Luo^,^ ^ 

^ Center for Interdisciplinary Studies & Key Laboratory for Magnetism and Magnetic Materials of the MoE, Lanzhou University, 
Lanzhou 730000, China 

^ Beijing Computational Science Research Center, Beijing 100084, China 
Received: date / Revised version: date 

Abstract. We investigate the pairing symmetry of the Kondo-Heisenberg model on triangular lattice, which 
is believed to capture the core competition of Kondo screening and local magnetic exchange interaction in 
heavy electron compounds. On the dominant background of the heavy fermion state, the introduction of 
the Heisenberg antiferromagnetic interaction {Jh) leads to superconducting pairing instability. Depending 
on the strength of the interactions, it is found that the pairing symmetry prefers an extended s-wave for 
small Jh and high conduction electron density but a chiral (1^2 _y 2 -h idxy-wave for large Jh and low con¬ 
duction electron density, which provides a phase diagram of pairing symmetry from the calculations of the 
ground-state energy. The transition between these two pairing symmetries is found to be first-order. Fur¬ 
thermore, we also analyze the phase diagram from the pairing strengths and find that the phase diagram 
obtained is qualitatively consistent with that based on the ground-state energy. In addition, we propose an 
effective single-band BCS Hamiltonian, which is able to describe the low-energy thermodynamic behaviors 
of the heavy fermion superconducting states. These results further deepen the understanding of the an¬ 
tiferromagnetic interaction which results in a geometric frustration for the model studied. Our work may 
provide a possible scenario to understand the pairing symmetry of the heavy fermion superconductivity, 
which is one of active issues in very recent years. 

PACS. PACS-71.10.Hf electron phase diagrams and phase transitions in model systems - PACS-71.27.+a 
heavy fermions 


1 Introduction 

The Kondo lattice model is at the heart position to un¬ 
derstand the ground-state property of the heavy fermion 
compounds determined by the competition between the 


the heavy fermion state of the Kondo-Heisenberg model 
on the square lattice is superconducting unstable and it 
favors a d-wave pairing symmetry, which is reminiscent of 
the d-wave symmetry of the Cu-based high-temperature 
superconductors resulted from the low-energy effective t-J 
Kondo effect and the Ruderman-Kittel-Kasuya-Yosida (RKK^^^del on the square Cu02 plane. [H] Due to the puzzling 


interaction (i.e., Doniach picture) both obtained in princi¬ 
ple from the intrasite Kondo coupling. [Tl[2l[3j|4] However, 
among the Ce-heavy fermion compounds, some exhibit 
antiferromagnetic long-range order, [5l[6l[71[8] which moti¬ 
vated Glesias, Lacroix, and Goqblin to revisit the “Do¬ 
niach picture” by including explicitly the intersite Heisen¬ 
berg interaction. pirniiTT] The revisited version of Doniach 
picture is able to capture the magnetic nature of those Ce- 
compounds. m 

The introduction of the Heisenberg coupling leads to 
many novel phenomena. naiiaiiaiiaiiTiiiH] Among those 
phenomena, the occurrence of the superconductivity is ap- 


pairing mechanism of the cuprates, the search of uncon¬ 
ventional superconductors becomes one of central topics 
in condensed matter physics in recent three decades. [22l 


On the other hand, the water-intercalated sodium cobal- 
tates Naa;Co 2 • ^H 20 was found in layered metal oxides 
Naa;Co 2 through a chemical oxidation process and the 
families of organic charge-transfer salts /t:-(ET) 2 X and 
Pd(dmit )2 are two examples to exhibit unconventional 
superconductivity. [77l[^[j^l30[l3 1 [[3^ The essential physics 
of these materials lies on the geometrically frustrated tri¬ 
angular lattice. j3^[3ll[35l[36[l37j The interplay of the elec- 


parently beyond the physics of the original Kondo lattice. [H tronic correlations and the geometric frustration is the 
[20] In a recent paper, m Liu, Zhang, and Yu found that physical origin of many exotic emergent phenomena like 

_ quantum spin liquid states. j38l[3^ In fact, it is also be- 

Present address: zhanglanl4@lzu.edu.cn lieved that the interplay could lead to a chiral pairing state 
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These recent developments motivate us to further explore where ct^(Qcr) denotes the creation (annihilation) opera- 
the superconductivity of the Kondo-Heisenberg model (KHM);or of the conduction electrons with spin cr. The first line 


on the triangular lattice. 

The triangular lattice has a Cqv rotational symmetry, 
and it allows a doubly degenerate E 2 representation of 
the superconducting order parameters with dj. 2 _y 2 and 


dxy degenerate states. 


This indicates that the KHM on 


the triangular lattice has a chiral superconducting ground- 
state and should exhibit a dx‘ 2 -y 2 + idxy{d + id) pairing 
symmetry. 

In the present work, we employ the large-N mean-field 
approach to study this model. As expected, we find that 
the superconducting state with chiral dEid pairing sym¬ 
metry can be possible for large Heisenberg antiferromag¬ 
netic interaction {Jr) and low conduction electron density 
{ric) from the calculations of the ground-state energy. Sur¬ 
prisingly, for small Jr and high Uc, the pairing symmetry 
is found to favor an extended s-wave. Thus we obtain a 
phase diagram of pairing symmetry of the model. Further¬ 
more, a qualitatively similar phase diagram has also been 
obtained from the calculations of the pairing strengths, 
which signals that the local Heisenberg antiferromagnetic 
interaction drives pairing instability with Fermi surface 
topology. This provides a physically intuitive understand¬ 
ing of the existence of the extended s-wave pairing sym¬ 
metry in such a heavy fermion system. It is also found that 
the phase transition between the extended s-wave and the 
chiral d+id-wave is a first-order when tuning Jr and Uc- In 
addition, we analyze the characteristics of the spectra and 
the density of states of quasiparticles for different pair¬ 
ing symmetries. In the superconducting state, an effective 
single-band BCS Hamiltonian is proposed to describe the 
superconducting properties such as superfiuid density of 
the model. Finally, one notes that the issue of the pairing 
symmetry of heavy fermion systems becomes very active 
due to some experimental and theoretical advances [42l|43l 
miii] and our work provides a possible scenario in un¬ 
derstanding the pairing symmetry of these heavy fermion 
systems. 

The paper is organized as follows. In Secl2j we intro¬ 
duce the Kondo-Heisenberg lattice model under the large- 
N mean-field theory. In Sec 131 we present the results and 
give the phase diagram of the pairing symmetries both 
from the calculations of the ground-state energy and the 
pairing strengths. The characteristics of the quasiparticle 
spectra, the density of states, and the pairing functions 
have also presented. Some discussion on an effective single¬ 
band BCS model is also given. Finally, SeclUis devoted to 
a brief conclusion and perspective. 


2 Model and mean-field approach 

The Hamiltonian of the KHM on the triangular lattice can 
be written as 


H = 






E ' 




-^E 

EJk Si • Si + Jr S 

* {ij) 


i(j 


( 1 ) 


in Eq. o describes the hoppings of the conduction elec¬ 
trons and /i is the chemical potential. The notations of 
(•) and ((•)) represent the nearest-neighbor (NN) and the 
next-nearest-neighbor (NNN) hopping respectively. The 
introduction of the NNN hopping term is to avoid the 
occasional nesting. The Jr (the coupling strength) term 
in the second line denotes the Kondo coupling between 
the localized f-electrons and conduction electrons. The 
Si = I fia^apfi/3 IS the Spin operator of localized 

f-electron with the local constraint — 1, the 

Si = ^ Scrcr' ^i opcrator of the con¬ 

duction electrons, and r is the Pauli matrix. The last 
Jr (the interaction strength) term is the Heisenberg ex¬ 
change interaction introduced by Coleman and Andrei 
to explore possible spin-liquid-stabilized unconventional 
pairing states. m It has also been used by Glesias et al 
to consider the antiferromagnetic long-range order in some 
Ce-based heavy fermion compounds. p fTQlITT] 

The standard large-N mean-field method [46] is suffi¬ 
cient to deal with the KHM on the triangular lattice. 
Since we consider the assumption that the pairing origi¬ 
nates from the antiferromagnetic exchange interaction, the 
spin-singlet pairing is expected to have lower energy com¬ 
pared to the spin-triplet pairing. In addition, the possible 
pairing symmetry should be compatible with the lattice 
symmetry of the triangular lattice. Thus, we consider two 
kinds of pairing order parameters, namely, the extended 
5-wave pairing Aij = Aq and the chiral d + id-wave par¬ 
ing Aij = breaking the parity and time-reversal 

symmetry, where Oij is the angle between the hopping di¬ 
rection and the horizon. The extended 5-wave has a uni¬ 
form phase Oij = 0, while the chiral d + id-wave has three 
possible values, namely, Oij = (0, |, ^). 

By the large-N mean-field, the Heisenberg antiferro¬ 
magnetic interaction can be decoupled by the particle- 
particle and particle-hole channels, i.e., Aij = ~{fi^f]i~ 

Vi = -(//,t/i,t + /k/i,i)> and the Kondo 
screening channel V = -h [H] As a re- 

suit, Si ■ Sj = - fljjA + iJ.c.] + = 


1 IXij (/if/it+/4/it) + 

field formulations, Eq. o can be rewritten in the k-space 
as follows 


t 


, 1^0' I 


. Based on these mean- 


k 


f Sk- 0 

0 —Sk + 

0 


V 0 


JkV 

2 


JkV 

2 

0 

Xk 


^k 




JhAI 
l ^JhAI 


0 \ 

JkV 

JhA}; 

-Xk J 


( 2 ) 


where = 


k - (4t’''-ki’4t’-^-k;) ^ four-component 


Nambu spinor, and Ak is the pairing function correspond¬ 
ing to the extended 5 -wave or the chiral d + id-wave pair¬ 
ing in momentum space. Sk = + 2ti [cos (^^/3ky) 
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+2 cos cos {^kx) with 7 /c = 2 cos cos (^) 

+ cos {kx) is the single-particle energy of the non-interacting 
conduction electrons. Xk = JnXlk + A is the kinetic en¬ 
ergy of the /-electrons, denoting the dispersion of real 
f-electron band. The parameter % = \Xij\ A is the 
Lagrange multiplier due to the local constraint of the lo¬ 
calized electrons. Ns is the number of sites in the lattice. 
The pairing order parameters for the extended s-wave and 
the chiral d + id-wdive read, 




cos{kx) + 2 cos 




= Ao 


cos{kx) 


COS 


cos 


cos 


2 






Diagonalizing the mean-field Hamiltonian, the quasiparti¬ 
cle energy spectra have two independent branches, which 
read 



Fig. 1. The ground-state energy Eg as a function of the 
Heisenberg antiferromagnetic interaction Jh in unit of t (here 
and hereafter the unit of energy is taken as t) for different 
conduction electron concentrations. The red dashed lines cor¬ 
respond to the extended s-wave pairing and the blue solid ones 
correspond to the chiral d + zd-wave pairing. The insets show 
the cases of small Jh- The parameters used are ti/t = 0.3, 
Jn/t = 2.5. 


where = U^k + + xl + , and £^2 = 

At zero temperature, the ground-state energy reads 


E„ 




^o^ , JkV^ , 3JhAI , 3JhX^ 
■ -E/c ) H ~ ^ ~ ^ ~ 


,_ (6) 

where E^ = E'^ + E^ = sj2 {Eki Ek2), ik = Sk - 

and Xk = Xk^k -^ 4 —• Minimizing the ground-state 

energy, one obtains a set of the self-consistent equations 
to determine the parameters x, Z\o, H, A, and the chemical 
potential ji: 


1 

1 

1 
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E 
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1 

Ek 

1 

Ek 

1 

Ek 

1 

Ek 


S/c 

Ek 2 






2 

Jk ’ 


Xk + 


^kXk 

Ek2 


= 0 , 


Xk + 


^kXk 

Ek2 


Ik = 3x, 


^ ^ Xk^kE^kJn^k 

Ns ^ Ek Ek2 

k 


1 - ric, 


(7) 

(8) 
(9) 

( 10 ) 

( 11 ) 


where the last equation is due to the constraint of the 
concentration of conduction electrons. 


3 Unconventional superconductivity 

3.1 pairing symmetry 

It is straightforward to solve self-consistently the mean- 
field equations given by Eqs. ©-(HU by assuming the 
corresponding pairing symmetry, namely, the chiral d-\-id 
(blue solid line) and the extended s-wave (red dashed 
line). First of all, we discuss the ground-state energies for 
these two different pairing order parameters. For different 
conduction electron density nc, the ground-state energy 
is plotted as a function of the Heisenberg antiferromag¬ 
netic interaction Jh, as shown in Fig. [TJ For all ric, it 
is found that irrespective of the pairing symmetry, the 
ground-state energy shows non-monotonic behavior of Jh, 
which indicates the competition between the Heisenberg 
antiferromagnetic interaction and the Kondo screening. 
Comparing with different pairing symmetry, it is found 
that for large Jh the system with the chiral d + id-wave 
pairing has a lower energy than that with the extended 5 - 
wave pairing. Surprisingly, with decreasing of Jh, the two 
lines go smoothly across, and as a result, the system with 
the extended s-wave pairing has a lower energy than that 
with the chiral d-hid-wave pairing. For very small Jh, the 
insets in Fig. [1] show that the ground-state energy with 
the extended s-wave pairing is still lower than that with 
the chiral d + id-wave pairing, although their energies are 
very close. This result indicates that for small Jh the sys¬ 
tem favors the extended s-wave pairing, rather than the 
commonly believed d-wave pairing in such heavy fermion 
systems. 

It is well known that the superconducting pairing sym¬ 
metry plays an important role in understanding of the 
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pairing mechanism. For example, for the weak-coupling 
superconductivity, e.g., the conventional BCS supercon¬ 
ductors, the pairing symmetry is an 5 -wave, in which elec¬ 
tron pairings due to the glue of phonon have a large weight 
on the on-site pairing amplitude in real space. For the 
strong-coupling superconductivity, e.g., high Tc cuprates, 
the situation is quite different. The strong Coulomb re¬ 
pulsion disfavors the on-site pairing but it drives the anti¬ 
ferromagnetic fluctuations. Thus, in high Tc cuprates the 
anisotropic pairing shows the d-wave symmetry. For the 
heavy fermion systems, many experimental and theoreti¬ 
cal results, e.g.. Ref. m, leaded to a commonly expecta¬ 
tion that the heavy fermion superconductors should have 
a d-wave pairing symmetry since the superconductivity 
is very close proximity to the antiferromagnetic region. 
However, some very recent works indicated that the 5- 
wave scenario is possible in the heavy fermion systems, 
which will be left for a brief review in the final section. 
Here we continue to discuss the characteristic behaviors 
of the heavy fermion superconducting states at different 
pairing symmetry. 

We first check the quasi-particle energy spectra ob¬ 
tained. One notes that nodal points will appear if Ek 2 = 0, 

which means that \Ak\ = 0 and need to sat¬ 

isfy at the same time. However, it is obvious that these 
two conditions are difficult to satisfy at the same time 
irrespective of the pairing symmetry assumptions of the 
extended s- or the chiral d -1- id-type pairing. 

It is found that the higher energy band branch of 
is away from the Fermi surface, which shows trivial fea¬ 
ture with different Jh- However, the energy band E^ is 
close to the Fermi surface. The numerical calculation of 
the quasi-particle energy spectrum E^ are shown in Fig. 
[ 2 ] for the d E id (blue solid line) and 5 -wave (red dashed 
line) pairing symmetry. They shows interesting features 
with different Jh- For small Jh = 0.075, the extended s- 
wave energy spectrum of E^ has two points very close to 
the Fermi surface which shows two negligible energy gap. 
In contrast, the chiral d + id-wave is almost gapless, as 
presented in the inset of Fig.[2](a). With increasing Jh up 
to Jh = 0.2292, the gap of the extended s-wave pairing 
increases, and a sizable gap begins to open for the chiral 
d+id-wave pairing, as shown in Fig.[2](b). This is the cross 
point of the ground-state energy for the extended s-wave 
and the chiral d + id-wave pairing symmetry at ric = 0.9, 
as shown in Fig. Hd). Further increasing Jh, both en¬ 
ergy gaps increase dramatically, and in this regime, the 
system always favors the chiral d + id-wave pairing sym¬ 
metry when comparing the ground-state energies of two 
kinds of pairing symmetry. From these discussion, one can 
conclude that the system favors an open but smaller en¬ 
ergy gap in heavy fermion superconducting state. A sim¬ 
ilar energy spectra feature has also been observed in the 
KHM on the square lattice. m Thus, it is reasonable that 
the extended 5 -wave pairing should also exist in the case 
of the square lattice. m 




Fig. 2. The quasi-particle energy band structure for the 
branch of E'J near the Fermi surface. The other branch of 
E'^ (not shown here) is away from the Fermi surface. The red 
dashed lines correspond to the extended s-wave pairing and 
the blue solid ones correspond to the chiral d -h zd-wave pair¬ 
ing. Here rzc = 0.9 and the other parameters used are the same 
as those in Fig. [T] 


as shown in Fig.[3l For small Jh, only the extended s-wave 
pairing shows a sizable energy gap, but for the chiral d-\-id- 
wave pairing assumption the energy gap is absent (Note: 
a finite density of states around the Fermi level for both 
cases is due to finite broadening in the calculations). This 
result is related to the above observation that the super¬ 
conducting states in the KHM always favors a finite energy 
gap, thus have the extended s-wave pairing for small Jh- 
With increasing Jh up to 0.2292 for Uc = 0.9[see. Fig. 
[3l[b)], the energy gap for the extended s-wave pairing as¬ 
sumption becomes more and more larger, and the chiral 
d + id-wave begins to open a sizable gap. After the quan¬ 
tum phase transition, for large Jh, the system begins to 
favor the chiral d + id-wave with a finite energy gap. A 
finite energy gap for the chiral d + id-wave pairing is con¬ 
sistent with the nature of the chiral d+id-wave pairing. [29] 
In addition, one also notes that in the chiral d + id-wave 
pairing case, the superconducting coherent peaks are ob¬ 
viously absent, which is sharp contrast to the case of the 
extended s-wave pairing. 

To further understand the superconducting instability, 
it is helpful to discuss the spinon pairing function in mo¬ 
mentum space, which can be deduced by using the stan¬ 
dard equation of motion method and the spectral theorem 
of the retarded Green’s function, as also done in Ref. m- 
The result is given by 



2^/2{Eki + Ek2) 



( 12 ) 


The transition of the pairing symmetry can also be re- The spinon pairing distribution for different Jh is shown 
fleeted by the conduction electron density of states calcu- by the first two columns in Fig. 0] for the chiral d + id- 
lated by using Green’s function G(/c, r) = — (T 7 -c/ccr(T)c^^( 0 )),wave pairing, in which the first column is the real part 



























Lan Zhang et ah: Extended s-wave pairing symmetry on the triangular lattice heavy fermion system 


5 



Fig. 3. The conduction electrons density of state for different 
Jh- The parameters used are the same as those in Fig. 


(with dx 2 _y 2 symmetry) and the second one is the imagi¬ 
nary part (with dxy symmetry), and the third column for 
the extended s-wave pairing. Comparison with these two 
cases, one notes that for small Jh, as shown in Fig. 0] 
(a),(d),(g), the pairing function of the chiral d + id-wave 
pairing is very weak but it is strong for the extended s- 
wave case. In this situation, the system prefers to be the 
extended s-wave pairing. With increasing Jh up to 0.2292, 
these two cases have comparable pairing functions, and an 
energy gap begins to open in the case of the chiral d + id- 
wave pairing. This is the point at which the system prefers 
to be the chiral d+id-wave pairing. Further increasing Jh, 
the energy gap of the chiral d + id-wave pairing continues 
to increase, but it is still smaller than the extended s-wave 
pairing case. Therefore, for large Jh, the system has the 
chiral d + id-wme pairing. 

We can also calculate the pairing functions for the con¬ 
duction electrons induced by the Kondo coupling. It is 
given by 


,t \ ^ JnJ^V^Ak 

8Ek2V^{Eki + Ek2y 


(13) 


As shown in FigISl crudely speaking, the pairing function 
of conduction electrons have the similar pairing symmetry 
as spinon, but show a 7r-phase shift in comparison to the 
case of spinon pairing shown in Fig 01 which is consistent 
with the result on square lattice. [16] It is noticed that the 
pairing intensity of conduction electrons is much smaller 
than its spinon counterpart, which reflects the fact that 
the pairing of conduction electron is induced by the pre¬ 
formed pairing of local electrons via non-vanishing Kondo 
screening. [ 1105 ] 

In the above discussion, we have fixed the concentra¬ 
tion of conduction electrons by Uc = 0.9. In order to study 
other concentrations of the conduction electrons, we have 
obtained a phase diagram of the pairing symmetry for the 
KHM model on the triangular lattice, as shown in Fig. 


Fig. 4. The spinon pairing distribution for different Jh in 
the first Brillouin zone marked by white lines. The first col¬ 
umn [(a)-(c)] is the real parts of the chiral d + id-wave pairing 
functions and the second column [(d)-(f)] is the corresponding 
imaginary parts. The third column [(g)-(i)] is the extend s- 
wave pairing functions. T, K and M are high symmetric points 
of the triangular lattice. The parameters used are the same as 
those in Fig. (T] 




Fig. 5. The conduction electrons pairing distribution for dif¬ 
ferent Jh in the first Brillouin zone marked by white lines. 
The first column [(a)-(c)] is the real parts of the chiral d + id- 
wave pairing functions, the second column [(d)-(f)] is the cor¬ 
responding imaginary parts. The third column [(g)-(i)] is the 
extend s-wave pairing functions. T, K and M are high sym¬ 
metric points of the triangular lattice. The parameters used 
are the same as those in Fig. [U 


[6][a) for Jk = 2.5. The phase diagram of pairing symme¬ 
try is separated by the cross points of the ground-state 
energy in Fig. (IJa). Due to the smooth cross of these two 
lines and different slopes, one can conclude that the phase 
transition is a first-order. Therefore, we do not expect rad¬ 
ical changes of thermal and transport behaviors near this 
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Fig. 6. (a) The phase diagram of the pairing symmetry for 

the KHM model on the triangular lattice in the ric — Jh plane 
for Jk — 2.5. (b) The phase boundary for different Jk- The 
parameters used are the same as those in Fig. [U 




Fig. 7. (a) The integrated pairing strengths of the extended 
s-wave (red dashed line) and the chiral d + id-wme (blue solid 
line) pairing symmetry and (b) the phase diagram of the pair¬ 
ing symmetry in the ric — Jh plane. The parameters used are 
the same as those in Fig. [T] 


transition point. For different Jk, the result is shown in 

Fig. mb). 

To understand the physics of the pairing symmetry 
transition, it is also useful to check the pairing strength, 
a quantity proposed by Hu and Ding in order to explain 
the resulting pairing symmetry on specific lattices. [l9] In 
the present case, the superconducting pairing strength is 
given by 


^ = ( 14 ) 

k 


where the form factor of the extended s-wave is /| = 
cos{kx) — cos (^) cos ^ ^ and that of the chiral d -f 

id-wave is = cos{kx) — iv^sin (^) sin — 

cos (^) cos ^ ^, The normal state quasi-particle ex¬ 
citation spectra are defined as ^ [{sk — + X/c) 


±\/(^/c - M - X/c)^ + {JrVY 


The result of the in¬ 


tegrated pairing strength is shown in Fig. [TJa). It is very 
clear that for small Jh and Uc = 0.9, the pairing strength 
of the extended s-wave is quite larger than that of the chi¬ 
ral d +id-wave. Increasing Jh, the pairing strengths show 
a cross point around Jh = 0.6355, which is identified as a 
first-order transition. When one changes ric, one obtains 
a phase diagram of the pairing symmetry according to 
the pairing strength, as shown in Fig. [TJb), which shows 
qualitatively the similar trend as the one obtained by the 
ground-state energy calculations, though the Heisenberg 



Fig. 8. The pairing strengths in the momentum space of 
extended s-wave pairing symmetry [(a) and (c)] and the chiral 
d + zd-wave [(b) and (d)] for different Heisenberg antiferromag¬ 
netic interactions. The black curve shows the Fermi surface. 
Here ric = 0.9 and the other parameters used are the same as 
those in Fig. [T] 


antiferromagnetic interactions have quite different values 
since the methods are different, as explained below. 

Fig. [5] shows the pairing strength in the momentum 
space and the corresponding Fermi surface at ric = 0.9 for 
different Heisenberg antiferromagnetic interactions. It is 
very clear that the pairing strength is quite sensitive to the 
shape of Fermi surface. On the contrary, the result from 
the ground-state energy is not so sensitive to the shape of 
the Fermi surface. Therefore, the phase diagram presented 
in Fig. 0 (b) can reflect more physics of the Fermi surface 
topology. In this sense. Fig. 0 (b) may be more practical 
for real systems, [29] for example, the water-intercalated 
sodium cobaltates Naa;Co 2 • ^H 20 . 


3.2 Effective single-band BCS description of 
superconducting states 

It is also interesting to note that although the KHM seems 
a two-band model, one may use an effective single-band 
BCS model to understand its basic features. [5Ql[5T| The 
main point is that for most of heavy fermion superconduc¬ 
tors, the superconducting state is believed to stem from a 
high temperature normal heavy Fermi liquid state, thus 
one can expect that the superconducting phase results 
from the pairing of renormalized heavy quasiparticles. 

Explicitly, we can first diagonalize Eq. m to obtain 
the heavy quasiparticles energy band using the Bogoliubov 
transformation C/(;cr — ^kPka, fka — ^/c^/ccr T^/c/^/ccr • 

Here the aj^a and Pka are the quasiparticle operators, and 
the Uk (vk) is the coherence factor. The Hamiltonian in 
terms of the new basis becomes 

H = E + EkJkJka) • (15) 

k(j 
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We have neglected the constant terms and set x = 0 since 
we only consider the pairing order parameter in Eq. m 
shown below. Then we add the pairing terms and project 
such terms into the renormalized quasiparticles basis to 
obtain the effective Hamiltonian as follows 


H = 


+Jh 

^Jh 


ka 


+ h.c. 
+ h.c. 


} 


(16) 


Since the Fermi surface lies in the /3-quasiparticle energy 
band, the pairing mainly appears in this band. Neglecting 
the (a-band pairing and the inter-band pairing, we can get 
an effective single-band BCS model for the heavy fermion 
superconductivity. 

H = [EkPl^Pka + [AkUkU-kPl^Piki + ^-c-) • 

ka 

(17) 

For this effective single-band model, the important trans¬ 
port measurement on the inverse squared magnetic pen¬ 
etration depth A will have ^ (x with ritot = + 1 

and being the effective mass of the heavy quasipar¬ 
ticle and Tic is the concentration of conduction electrons. 
The superfluid density experiment would observe the ex¬ 
pected exponent behavior at low temperature compared 
to the superconducting transition temperature since both 
the extended s-wave and the chiral id-wave are gapped 
and have no nodes for the present triangular lattice. In 
addition, one may also expect the optical conductance is 
gapped due to the pairing gap. However, it should be em¬ 
phasized that both magnetic penetration depth and op¬ 
tical conductance experiments are not able to distinguish 
the pairing symmetry between the extended s-wave and 
the chiral d + id-wave. Therefore, a phase sensitive ex¬ 
periment to confirm the true pairing symmetry is needed 
if some real-life heavy fermion materials are indeed de¬ 
scribed by the present models with triangular lattice. 


4 Conclusion and perspective 

In summary, we have investigated the KHM on triangular 
lattice with the fermionic large-N mean-field theory. A 
qualitatively same phase diagram of the pairing symmetry 
has been obtained by comparing the ground-state energy 
and analyzing the pairing strength. At small Heisenberg 
antiferromagnetic interaction and large concentration of 
conduction electrons, the system prefers to an extended 
s-wave pairing symmetry and on the contrary, the system 
favors the chiral d + id-wave pairing symmetry. The phase 
transition between these two pairing symmetries is found 
to be first-order. 

In very recent years, the issue of the pairing symmetry 
of the heavy fermion systems becomes quite active due 


to some experimental and theoretical advances. For ex¬ 
ample, the specific-heat measurements in CeCu 2 Si 2 along 
with its linear dependence as a function of magnetic field 
and the absence of oscillations in the field angle suggest 
that this material may have an s-wave pairing symme¬ 
try. m The s-wave pairing symmetry has also been pre¬ 
dicted by a first-principle calculation in CeCu 2 Si 2 . [52] For 
another heavy fermion compound CeCoIns, the London 
penetration depth measurements also suggest that a T6 
substitution can change the pairing symmetry from nodal 
to nodeless, which is closely related to the Fermi surface 
topology change. [JS] When it is still debated what pairing 
symmetry the heavy fermion superconductivity is, j^l44l 
[45] our work provides a possible scenario to understand 
the pairing symmetry in heavy fermion superconductors. 
To test our result, it is hoped that future experiments 
can discover some heavy fermion compounds with the ge¬ 
ometrical frustrated triangular lattice structure, which is 
highly desirable from the interest in unconventional super¬ 
conductivity and the general frustration/interaction heavy 
fermion phase diagram. [53] 
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